An Extension of the Van der Waals 
Equation of State 



Jianxiang Tian 

Department of Physics, Dalian University of Technology 
Dalian 116024, P.R.China 
Department of Physics, Qufu 
Normal University, Qufu 273165, P.R.China 

Yuanxing Gui 

Department of Physics, Dalian University of Technology 
Dalian 116024, P.R.China 

February 2, 2008 



Abstract 

In this paper, we extend the Van der Waals equation of state to a 
universal form and find that the Van der Waals equation of state is a 
special condition of this form. Then a good form of equation of state for a 
balanced liquid-gas coexistence canonical argon-like system is given. The 
localization of this universal form is mentioned out in the end. 
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1 INTRODUCTION 

Much attention has been paid in recent years to the hard core Yukawa (HCY) 
potential as a model for the pair interactions of fluids [1]. The liquid state theo- 
ries such as the Mean Spherical Approximation (MSA) and the Self Consistent 
Ornstcin-Zcrnikc Approximation (SCOZA)are proposed. Recent studies of the 
HCY fluid can be found in [2,3] and references therein. But the Van der Waals 
equation of state 

v — b v 

was less mentioned recently. Here a and b are parameters, ks is the Boltzmann 
constant, and v — V/N. In PART 2 of this paper, we extend the VDW equation 
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of state to a universal form and find that the VDW equation of state is a special 
case of this form. In PART 3, we state the localization of this universal form 
by applying it to the liquid-gas phase transition to one order. 



2 Theory 

We define particle number density by 



_ N _ 1 
V v ' 



(2) 



Then Eq.(l) can be expressed as 



n Nk B T 2 
P= V-Nb- an - (3) 

Now we extend the VDW equation of state to such a form as 

B and a are parameters, too. When a = 1, B = a, Eq.(4) reads the VDW 
equation of state. Now we sign 

T* = T/T c , (5) 

n* = n/n c , (6) 

P* - P/Pc, (7) 

with T c the critical temperature , P c the critical pressure , n c the critical particle 
number density. 

At critical point, the function P = P(V,T 7 N) has such qualities as 

dP 

^k=0, (8) 
d 2 P 

Thus we get the critical data from Eq.(4) by solving Eq.(8) and Eq.(9). They 
are 

k B T c = <j(a + l)(^—) 2 Bn° c , (11) 
a + 1 

P c =^-Bn°+\ (12) 
a + z 
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_ n c k B T c _ 4{a + 1) 

(13) 

Letter C represents the critical coefficient. Substituting Eq.(10-12) to Eq.(4), 
we get the reduced equation of state: 

p . = 4n*T* (a + 1) {a + 2)n*^ 

((ct + 2) -n*a)a a ' 

When (j=l, Eq.(4) reads 

which is just the form of the VDW equation of state. Data in Eq. (10-13) read 





(16) 




(17) 


M = 276' 


p - s 

c 276 2 ' 


(18) 


n c k B T c 

p ^v 3 - 


(19) 



We are very familiar with these results, which are the data from the VDW 
equation of state. Eq.(14) reads 

<S " 7 " in* 2 , (20) 



(3-n*) 



which is the reduced form of the VDW equation of state . From Eq.(3), Eq.(4) 
Eq. (14) and Eq.(20), we see that the original form of the VDW equation of 
state is 

Nk B T 1+1 
P= V-Nb- 1 * an ■ 



(21) 



When cr = 0.7432, Eq.(4) reads 



P = 1]^L _ .7432Sn 17432 . 
V - Nb 



(22) 



Eq.(10-13) reads 



n c = 0.27096 -1 , (23) 

k B T c = 0.2609S6- 07432 , (24) 

P c = 0.0207S6- 17432 , (25) 

C= --f,3.4201, (26) 
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Eq.(14) reads 



9.3821n*T* t1 , N 

-3.6911n* 17432 . (27) 



(2.7432 - 0.7432n*) 



Eq.(26) offers the critical coefficient of argon approximatively. 

In 1975, W.G. Hoover, etc, generalized the Van der Waals mean-feild attrac- 
tion by choosing the attractive potential energy proportional to the Sth power 
of the density [4]. Its per-particle partition function has the form[4]: 

V " cxp [_%__M_ + _^_^]. (28 ) 



NX" Nk B T Nk B T k B T 

Parametres are interpretated there. S is the exponent of n. We all know that 
the equation of state can be obtained by 



Then we can integrate Eq.(4) to get its corresponsive per-particle partition 
function. The result is 

z = f(T,N,m)(V-Nbfexp(^ f n^ , (30) 

with / a function indenpendent with V. Here m is is the mass of one particle. 
The equation of state obtained from Eq.(28) is complex. It is difficult to get its 
reduced form. Relatively, Eq.(4) is easier to be operated. 

If we replace the Van der Waals repulsion term in Eq.(4) with the 
Carnahan-Starling hard-sphere repulsion term[5,6], the eqution of state is writ- 
ten as 

Ay - 2y 2 



Nk B T 
V 



1 + 



-aBn a+1 , (31) 



with y = gj {AV m ). Here g is the covolume of Carnahan-Starling hard-spheres [5], 
and V m is mol— volume. The critical data of Eq.(31) and the reduced form of 
Eq.(31) are difficult to be solved. We do not discuss it in this paper. 



3 Discussion 



Here we only study the balanced liquid-gas phase coexistence canonical system 
with Eq.(4) and Eq.(14). 

In 1945, E.A.Guggenheim collected the data of this system from experiments 
and gave out the correlation of p*, p\ and T* by the empirical equations [7] below 

p* = l + 0.75(1 - T*) - 1.75(1 - T*) 1/3 , (32) 
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For 



p* 2 = 1 + 0.75(1 - T*) + 1.75(1 - T*) 1/3 . 



(33) 



p\ = mni/mn c — n\/n c — n*, (34) 

p2 = mn 2 /mn c = n 2 /n c = n 2 , (35) 

we have 

n* = 1 + 0.75(1 - T*) - 1.75(1 - T*) 1/3 , (36) 

n* = 1 + 0.75(1 - T*) + 1.75(1 - T*) 1 / 3 . (37) 

In Eq. (32-35), p\ is the reduced density of the gases and p 2 is the reduced density 
of the liquids. As far as an argon system is concerned, the inaccuracy of these 
two equation is generally only one or two parts per thousand of p 2 or of p\ when 
T* > 0.60T C [7]. So, it is acceptable to consider the data of (T* ,nl(T*),n* 2 (T*)) 
from Eq.(36) and Eq.(37) as the experimental ones in this temperature region 
when the balanced liquid-gas coexistence argon-like system is considered. In 
Eq. (34-35), m is the mass of one particle. E.A.Guggenheim gave a numerical 
analytic result of the relation between the reduced temperature and the reduced 
pressure from experiments by equation[8] 

P* = exp(5.29 - 5.31/T*), (38) 

which best fits the experimental data for argon when T* > 0.56T C except a 
tiny region near the critical point [8]. Thus it is acceptable to consider the data 
of (P*,T*) from Eq.(38) as the experimental ones, too. Substituting the data 
(T* ,n*(T*)) from Eq.(36) to Eq.(14), we get the data of the reduced pressure 
of the gases when a is fixed. Then we compare these data with the ones gotten 
from Eq.(34) to see the accuracy of Eq.(14) when a is chosen to be this value. 
Table. (1) illustrates the data with a = 0.7432, a = 0.8, a = 1, and a = 1.5. Pj" 
is the reduced pressure of the gases gotten by substituting the data (T*,n*(T*)) 
in Eq.(36) to Eq.(14). And P 2 * i s the reduced pressure of the liquids gotten by 
substituting the data {T*,n*{T*)) in Eq.(37) to Eq.(14). From Table 1, we can 
see 

• |pr-p;u=o.74 32 < |pr-p e *| CT =o. 8 < |p;-p;l=i < |pr-p;| CT= i. 5 when 

T* > 0.8. 

• I-P2* ~-P e *U=o.7432 < |P 2 * — P*|cr=o.8 < |P 2 * -P e *|(T=i < |P 2 * -P*|(t=i.5 when 
T* > 0.8. 

• |Pj* - P*| ctT * < |P 2 * - Pe\aT> when T* are fixed in the region 0.7 < T* < 
0.99. 

So we can conclude that a — 0.7432 is the best one of these four values in 
the temperature region > T* > 0.8 and Eq.(14) fits the gas phase better than 
the liquid phase. 
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Here we will ask whether a proper a = <jo, which can bring right academic 
forecast to the experimental data P* , n*, n 2 and T* , exists. Now we suppose 
that such a proper cr exists. Then Eq.(14) is fixed to be the form of 



p* 



An*T* (a + 1) 



(go + 2) n*( g " +1 ) 



(39) 



((cr + 2) - n*er )o-o 



Thus we have 



2/ = 2/(J ,ct ) = P ,f7 ) - P (n 2 , T >^o) = 0, 



(40) 



where n*, n 2 arL d ^* are the experimental ones in Eq.(36) and Eq.(37). Accord- 
ing to our supposion, at two different reduced temperature Tj" and T 2 , Eq.(40) 
will be true. Fig.(l) is the curve of y = y(T£ = 0.8, a). Fig. (2) is the curve 
of y = y(T 2 * = 0.85, a). From the figures, we know that y = y(T£ = 0.8, 
a) = when a = oq(T{ = 0.8) = 1.7050, y = y(T 2 * = 0.85, a) = when 
a = cr (T 2 * = 0.85) = 1.9190. For er (?I = 0.8) ^ (T (T 2 * = 0.85), the proper 
a = (Jo we try to find does not exist. 

In 1990, J.M.Kincaid, etc, representated the pure fluid coexistence curves 
by series expansion successfully [9]. Detailed description about pure fluid coex- 
istence curves such as argon's can be seen there. Our result here is not as good 
as theirs, but is more simple in calculation. 
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Figure 1: The curve of y = y(T£ = 0.8, a) versus a. The dotted line is beeline 
y = which is a referrence here. 
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Figure 2: The curve of y = y(T* = 0.85, a) versus a. The dotted line is beeline 
y = which is a referrence here. 
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